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Abstract 

We discuss certain types of cyclic and nearly cyclic interactions among ./V "point" - 
agents in the plane, leading to formations of interesting limiting geometric configura- 
tions. Cyclic pursuit and local averaging interactions have been analyzed in the context 
of multi- agent gathering. In this paper, we consider some nearly cyclic interactions that 
break symmetry leading to factor circulants rather than circulant interaction matrices. 

1 Introduction 

Consider a "swarm" or "pack" of N robots in the plane, denoted by Vo,Vi, ■ ■ - Vn-i 
which can all see each other and are aware of the other robot's identities (i.e., can 
distinguish them). We shall define the rules of interaction specifying how each robot Vk 
moves in response to the (evolution in time of the) configuration of the entire swarm. 
Therefore denoting TVs location at time t to be Vk(t) = Xk{t) + iyk{t) (a complex 
number), we assume that we can write the swarm evolution equations as follows: 

or V k {t + l) = $l D) {V s (0\s=o,l,...N-i;£<t} (1) 

depending on whether the temporal evolution is continuous (C) or discrete (D). So far 
the ^-operators are not specified, and in fact they could be quite involved in general. 
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(C) 

The operator & k provides an instantaneous velocity vector for agent Vk m response to 

the locations of the other agents in the swarm, while $!l will yield the next location for 
Vk in a synchronous discrete timed evolution. These operators should produce the same 
motion if we decide to look at the agents in different frames of reference, i.e., re-encode 
their locations using transformed coordinates, hence the resulting equations should be 
at least similarity invariant, and maybe even affine invariant. The requirement to have 
the same evolution equation in arbitrarily similarity (i.e., scaled Euclidean) or affine 
transformed coordinates clearly imposes restrictions on the # operators and some of 
these will be discussed in the sequel. 

An important class of operators are the linear memoryless ones which have the form 

N-i 



$k{Po,Vi,...V N -i} = J2 m i( t )' p i( t ) 



1=0 

where m k {t) are some (complex) numbers, varying perhaps in time. In this case, 
Equation (pQ) describes a linear (generally time varying) system's state evolution, and 
there is a wealth of theory dealing with such systems in the control and signal processing 
literature. Here we shall mainly be concerned with a special class of (constant) linear 
Toeplitz operators of the form 



$*{Po,Pi,...Pjv-i} 
where A is some complex number, and 



jV-1 

£ A Ind[G-fc)<0] 
1=0 



m (l~k)mod N^lit) 



(2) 



m_i = mjv-imod iV 
m_k = m-jv-fcinod N 



and Ind[a < 0] 



1 if a < 
if a > ' 



Writing out explicitly Qki'Po, ■ ■ ■ , 'Pn-i} f° r k = 0, . . . , N — 1 in matrix form and 
denoting 



P(i) 



the swarm's evolution dynamics becomes 



|p(t) or)P(/ + i ) 



$P(t) 



(3) 



m 


mi 


m 2 


mAr_i 




m 


mi 


m N -2 


\m N _ 2 


Amiv_i 


m 




Ami 


Am 2 




XrriN-i mo 



P(t). 



Note here that if A = 1, the matrix is a special Toeplitz-circulant matrix, otherwise it 
is a generalization of a circulant called a A-factor, or A-circulant matrix. Such matrices 
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arise in several applications, such as linear systems theory 0[1O], linear algebra [T], 
geometry [5j [131 114] . and in connection with inverses of Toeplitz matrices [91 H], 
coding theory [6] and linear systems of differential equations [T7J. In case of A = 1, i.e., 
when the operator $ is Toeplitz-circulant, we have that all the robotic agents perform 
"cyclically" the same operation, i.e. agent Vk will determine its next location (or its ve- 
locity) according to the same weighted average performed on Vk-,Vk+i-, ■ ■ ■ 'P(k+N)modN 
(in this order), i.e. 



V k (t + 1) 
or iV k (t) 



which can be rewritten as 



[m ,mi, . ..m N . 



V k (t) 
V k +i(t) 

_ 'P(k+N)mod AT CO 



(4) 



V k (t + 1) 

A' 
dt' 



or iV k (t) 



m 



1 ( kth P lace ) ... o 



1 



mZ 





fc-i 



where 



1 

: 



1 



■P(t) 
" 



1 

1 




(5) 



and m = [mo, mi, . . . m^-i}. 



This special case, with a circulant matrix $, was extensively analyzed before in the 
context of polygon smoothing evolutions and cyclic pursuits for robotic gathering and 
formation control, see e.g. 0[I3O[3m[l2l[IH[7]. 

Note that invariance requirements impose some conditions on the linear evolution 
operators, as we now discuss. If P(t) is described by the evolution equations 



or P(i + 1) 



$(C)p(i) 
$(°)p(t) 



from some initial location P(0) = P(i = 0), and if we re-encode the agents' positions 
via a general similarity transformation of the form 

where p and r are some complex numbers and 1 = [1, . . . , 1] T , we shall have for P'(t): 
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in the continuous case 

|p<( f ) 4 |(,P(*) + rl) 

which is equal to (pP(t) + rl) only if & c h = 0. 
in the discrete case 

P'(t+1) = pP(t + l)+rl 
= p$( D )p(t)+rl 



which is equal to $( D ) (pP(t) + rl) only if ¥ D h = 1. 

Hence the ^-matrices that describe linear, time- invariant evolutions need to obey 
the conditions $(^1 = or = 1 in order to have Euclidean or similarity 

invariant evolutions. In some of our examples, these conditions cannot be satisfied. 
However, note that any N x N matrix $ may be embedded in an (TV + 1) x (N + 1) 
matrix 3> as follows 



$ s 




' 1 " 




" $1 + s " 


z 




1 




z 



and selecting either z = and s = — <3?1 or z = 1, we obtain a $ matrix that describes 
an invariant evolution of a multi- agent system with an additional agent Vb whose 
position is stationary (^Vb = or Vs(t + 1) = ^s(t)). This additional agent will act 
as a "beacon" or a set reference point, for the description of the swarm of agents. In 
this case, setting Vb = (0, 0), the evolution of the rest of the agents will be described 
by the original matrix Note that the spatial location of the fixed Vb in the plane 
may be determined according to the initial location of the agents of the swarm. A good 
example is the geometric and affine invariant decision that can be made by each agent 
independently to set Vb, and hence the origin of its Cartesian coordinate system, at 
the centroid of the agent location constellation at t = 0. This will make the swarm 
evolution entirely autonomous. However, an external setting of the location of Vb 
might be useful in controlling the swarm and steering it toward a desired place in the 
environment. One might even desire to move Vb in time and make the swarm move 
accordingly, by tracking the beacon point in addition to its own internal dynamics 
controlled by <I>. 



2 Analyzing Swarm Evolution via Mode Decou- 
pling 

Circulant, and A-factor circulant matrices have very special structures and this allows 
us to diagonalize them, essentially by Fourier transform methods. Let us see, in general, 
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how diagonalization yields a way to analyze the evolution of the constellation of robots 
by decoupling it into independently evolving modes. Indeed assume that the time- 
invariant matrix $ can be diagonalized (for example when has distinct eigenvalues, 
hence a full set of orthonormal eigenvectors), as follows 

$ = T~ l DT 



where D = Diagfdo, d± . . . djv-i] displays the eigenvalues of $ and the columns of T 1 
are the (right) eigenvectors. Now we have that 



or 



P(t + 1) 



T^DTVit) 



and hence 



TP(t + 1) 



D(TP(t)). 



In terms of the transformed vector P(t) = TP(t), the evolution is a decoupled evolution 
controlled explicitly by the (constant) eigenvalues [10J. Indeed, we have 



P(t) 



d\ 



a N-l 



P(0) 



(discrete case) 



or 



P(t) 







P(0) 



(continous case) 

Therefore diagonalization enables the explicit solution of the swarm evolution, in the 
case the matrix is time invariant and has a full set of orthonormal eigenvectors. As 
we shall see below, A-factor circulants are a family of matrices that enable both a nice 
physical interpretation in terms of cyclic and symmetric interactions among similar 
agents and an explicit diagonalization via discrete Fourier transform matrices. 



3 Diagonalization of Factor Circulants 

Factor circulant matrices are very special in that they provide explicit formulae for 
the diagonalizing transforms and for their eigenvalues. This enables us to analyze in 
detail the behavior of multiagent interactions when these are cyclic or "nearly" cyclic, 
and fully describe the limiting behaviors of the swarm. For circulants, we have the 
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following results. Consider the unitary Fourier transform matrix 

a ( ; ... ( ' ' ' 

n i 

1 



FT 



N-l 



w° w N 1 



(JV-1)(JV-1) 



1 



in 



(fe-i)G-i) 



k,l=l,...,N 



where w 
only if 



e 1 n is an N root of unity. Then C is a Toeplitz-circulant matrix if and 

Mo 



CfFTl = fFTl 



Mi 



/xjv-i 

where /Uoi Mi) • • • > MiV-i are the eigenvalues of C and are given by 

N-l 

w = 



k=0 



Hence 
and 



[FT]*C[FT] = Diag[^o,Mi,--- ) MJV- 



C = [FT]Diag[Mo, Mi, ■ • • , MJV-i] [FT]*. 
To summarize the remarkable properties of circulants, we can state that they are (1) 
diagonalized by the discrete Fourier Transform, (2) they all commute, (3) their products 
are circulants, (4) their sums are circulants too, and (5) their inverses/pseudoinverses 
are circulants, and are readily found [9]. In fact, many of the wonders of modern signal 
processing algorithms, and linear, time invariant systems theory stem from the above 
properties. 

The corresponding, and equally remarkable properties of A-circulants are, however, 
much less known and applied. Suppose we consider the following operation on a circu- 
lant C = Cr, 



|co,ci,...,cjv-iJ ' 



W 



di 







ajv-i 



[c ,ci,...,cjv-iJ 







b N . 



i.e. W is obtained by pre- and post multiplying C by two diagonal matrices. It is easy 
to see that we have 



W= C 



C0,C1,...,CJV- 



.1]© 



a b 
aibo 



aobi 
a\bi 



■ ■ ■ aobjsi-i 
. . . aibN-i 



COM 
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where stands for the Schur Hadamard multiplication (or a "masking" operation) 
which multiplies matrices element-wise, and 

M = [a k bi] k ,i=o,...,n-i- 

For matrices of the type W, we have that they inherit interesting diagonalization 
properties from the original circulant C. The matrix W is a circulant matrix that is 
modified by a highly structured masking matrix M and we have that 

W = Diag[a , . . . , ajv-i] [FT]Diag[/i , . . . , /ijv_i][FT]*Diag[6 , . . . , 6jv_i]. 

However, since the masking matrix is neither circulant nor Toeplitz, we shall have to 
consider some special cases for the {ao, ai, • • • , ajv-i} and {60, bi, . . . , bN-i} sequences. 
First of all, note that the factorization above will be of the form 



W = U 



Mo 



if and only if 



(Diag[a ,ai,...,ajv_i][FT]) 1 
W-Diag/;,. 1 - 



I^N—l 



u 



[FT]*Diag[6o,6i,...,6iv-i] 
[FT]*Diag[6o,6i,...,6iv-i] 



or bk = a k 1 , and U will further be unitary if also bk = a* k , implying that a& = e^ ak and 
bk = e _jafc = a* k . In this case the masking-matrix multiplying C will be [e^ ak e~^ ai ] = 

[etf a *- Q ')] fc ,,= 0) ...,jv-i- 

The most interesting particular cases of {ao, a±, . . . , oat-i} and {bo, bi, ... , 6jv-i} 



arise when we have 


ak = 


7 fc and b 


k = 7 


k , k = 


0,1. 


..,N- 


-1, for 


some 


7. In this case, we 


have 


in general 
















" 1 


7- 1 


7" 2 




7 -(JV- 


1) " 








7 


1 


7- 1 










M 




7 2 


7 


1 




7 -(JV- 


l)+2 








yv-i 






7 


1 
















1 


1 


1 


1 1 












7^ 


1 


1 


1 1 




= Circ [1>7 _i 


...,7-(w 


-i)] 


1 N 


7^ 


1 


1 1 



7 ^ 7 ^ 



where Circ^-i^ __ j7 _ (JV _i)] is given by 



1 

7 -(jv-i)+i 7 



-1 



-2 



7 7 

1 7 -i 
-(;v-i) j 



7 



7 ^v ! 



7 -(JV-i) 

-(JV-1)+1 



7 



-(JV-l)+(JV-2) 



1 

-(AT-1) 



7 



Hence the matrix W = C M becomes 



W = C [c0j ... jC ^_ l] Ckc^-!^-^-!)} 



11111 

7 "1 1 11 
yv yv 1 1 1 



<y N 7^ ... 7 JV 1 



N 



which clearly is a A(= 7^)- circulant matrix. 

To summarize, we have the following result: A A-circulant matrix W, denoted by 



W 



mo mi mi 
Xmjy-i mo m± 
Amjv-2 Amjv_i mo 



m N -i 
m N - 2 



Ami Arri2 • • • AmAr_i mo 
can be rewritten as 

W = Circ [m0)mi7jm272j ... )mjv _ l7 jv-i] Circ [li7 _i 5 ^ )7 - ( jv-i)] A 

with 



A 



and hence can be factorized as 



1 


1 


1 




1 


A 


1 


1 




1 


A 


A 


1 




1 






A 


1 




A 


A 




A 


1 



and 7^ = A 



W 



7 



JV-l 



[FTl 



Ho 



Hi 




H2 



Hn-i 



[FTl 



7 



7 



-(N-l) 



where [/xq, Hii • • • , /-tjv-i] are the eigenvalues of 



given by 



Ci rc [mo,mi7,...,mj V _i7 JV - 1 ] ^ tc [co,ci,...,cn-i] 



N-l 



Hi = ^m fc -7 fe -e (7 = A ^ ) . 

fc=0 

Therefore W is readily diagonalized as follows 



1 



1 



Ho 



7 



-1 



7 



[FTl 



7 



W 



7 



[FTl 



Hn-i 



-(N-l) 



7 



7V-1 



= T X WT, 
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the matrices T and T 1 being 
1 

7 



7 



JV-l 



[FT] and T" 1 = [FT] 



7 



7 



JV-l 



Note that T is not, in general a unitary transformation. In all developments above, 
we assumed 7 to be arbitrary. If 7 7^ is a real number, T will be an invertible 
matrix, as seen before. If however 7 is purely imaginary, i.e. 7 = e Jlfi , then clearly 
7* = e~^ v = 7 _1 and the matrix T becomes a unitary transformation, obeying 



T*T 



In this case the matrix W will be A-factor circulant with A = e^ N . 



4 Dynamics of a Cyclically Interacting Swarm 

Returning to the problem of analyzing the dynamics and the long-term behavior of a 
swarm of robots Vq,Vi, . . . ,Vn-i interacting according to 



P(t + 1) 
or -P(t) 



mo mi 777,2 

AmjV_l 7770 "71 

XrriN-2 Amjv-i mo 



m N -i 



P(t) 



A?77l A7772 . . . Am/V-l 7770 

= *P(*)> 

we have that the interaction matrix <I> is A-circulant hence it is diagonalizable as follows: 



7 



7 



JV-l 



[FTl 



Ho 



Hn-i 



[FTl 



7 



7 



7" 



-(JV-l) 



where 7 = A^ and 



Therefore defining 



JV-l 
k=0 



P(t) 4 [FT]* 



A jv 



A jv 



P(t) 
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we have decoupled dynamics for the transformed location vector, given by 



|P(*) 
or 

P(t+1) 



/'o 



Ml 




fJ'N—l 



P(t) 



and the evolution of the swarm is controlled by the eigenvalues hq, fj,±, . . . , //jv-i- 

Let us concentrate next on some specific cases of m = [mo, . . . , mjv-i] and A. A 
"A- cyclic" interaction involves agents that are reacting differently with the agents that 
follow them to the agents that precede them in the ordering Vq, ■ ■ ■ ,Vn—1- 



4.1 Darboux's polygon evolution and extensions 

As a first example, suppose that we have a generalization of Darboux's polygon evolu- 
tion process [5], which is also a nice model for cyclic pursuit: 

5 1° ... 
1 



P(t + 1) 








A± 



1 l 

2 2 



P(t). 



In this case, we have a A-factor circulant with 

1 



1 1 ,i _,-2zl.; 

- + -\Ne N 



-(1 + Aive 



Z = 0,l,...,iV-l. 



Here, the evolution of the polygon vertices (or the agents in cyclic pursuit) is described 
by 



P(t + l) 



where we defined 



P(i) 



fFTl 



Mo 

Mi 






Mjv-i 



P(0) 



^-(JV-l)/JV 



P(t). 
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From this we have 



P(f) 



A 



[FT]P(t) 



N-l 

A— 



[FTl 



Mo 



o '•• 



I 1 N-l 



P(0). 



The evolution of the polygon vertices (the swarm of robots) when we let the time 
grow, thus asymptotically depends on the dominant eigenvalues among fiQ, . . . , //jv-i- 
In the case of A = 1 (or circulant cyclic pursuit), we have 



w = I(l + e -^), Z = 0,1,...,7V-1, 



and jUo = 1. Then 

P(*)t-»oo 



[FT] 



m 



[FTl 



Mi 
'■. 

Mi 



Mjv-i 



P(0) 



t^oo 



Mjv-i 



[FT]*P(0). 



t— »oo 



Since the dominant eigenvalue /Uo = 1 and all others have modulus less than one, we 
have that the limiting behavior is 



P(t)t->oc = ^ 



[1,1,...,1]P(0). 



Hence the point constellation converges to the centroid of the initial locations. The 
way this convergence occurs is be controlled by the next dominant eigenvalues, which 
are in this case 

1 -,'2ZLx 

Mi = ^ 1 + e ^ 

1 . 2tt(JV-1) 

Miv-i = 2^ 1 + e 1 N ^ 
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Indeed, writing 



P*(i) = P(t) - i 



[1,1,...,1]P(0), 



we have 



P N (t) = [FT] 



Ml 

'•• 



Miv-i 



[FT]*P(0) 



and, disregarding the faster decaying terms /4, « = 2, . . . , N — 2, we further get 



P*(t)^oo 



1 

N 



w 



1 

w 

N-l 



[i,w,...,w N - 1 ]-p(oH 



+ 



N 



1 

N-l 



w 



w 



(N-1)(N-1) 



Hence 



1 

N 



w 



1 

w 

N-l 



w 



1 

AT-1 



W {N-X){N-1) 



£(*)Mat-i 



where A(t)/j,\ and .B(i)/ijv-i are some complex numbers, and P N (t) will be, in the 
limit t —* oo, an affine transformation of a regular polygon, i.e. a discrete ellipse (see 
Figured]). 

For the general case where A is some real or complex number, we have that 



P(t) 



t^oo 



N-l 

A jv 



[FT1 



Mo 



JV-1 

A 



[FT1 



Mo 



Mi 






o '•. 



P(0) 

P(0), 

t— >oo 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.6204 0.6204 0.6204 0.6204 0.6204 0.6204 0.6204 



Figure 1: The cyclic pursuit case (A = 1) with a random initial polygon with 
N = 7 points, the first figure presents the initial configuration (in red) and the first 
iteration (in blue), the second shows the entire evolution for 100 iterations, the last 
figure displays the scaled up configuration for the last few iterations. 



where [Iq = h(l + A 1/W ) is the dominant eigenvalue. Since 



P(0) = [FT] 



X-l/N 



-(iV-l) 
A N 



P(0), 



we then have that 
1 

P(£)t->oo = fJ-o 



N-l 
A N 



[FT] U ([FT] 



1,1, 



A-i/JV 



-(JV-l) 
A N 



P(0) 



and since the first column of the Fourier transform is a vector of all ones, this further 
simplifies to 



P(t) 



t^oo — Ho 



N 



1 



JV-l 
A N 



1 IAT -(JV-l) 

^A-V^^.^AT^O). 



Therefore, we see that the limiting behavior is dominated by 



P(t) 



t— >oo 



1(1 + aH 1 [i,a-^,... , a^]p(o) 

2 iv > „ ' 



a (complex) scalar 
We can distinguish different behaviors depending on A. 



1 



A n 
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Figure 3: A = — 1 



1. if A is real and |A| < 1, P(t) tends to zero, but the limit behavior will be a linear 
constellation of points 



" 1 " 




' 1 " 


A l/iV 








+ i{a t ) y 










A N 




\ N 



If |A| > 1, the constellation of agent locations will diverge in a similar formation. 

2. If A is a complex number P(A)e l</3(A) , the convergence/divergence will depend on 
the angle of rotation induced by tpr^) and on the magnitude /3n\ . As seen in the 
examples provided in Figures El [3l |H [5l [H [7J in the limit, agents are marching in 
elliptic or circular arcs, spiralling towards their point of convergence (and in case 
of divergence, spiralling out to infinity). As in Figure [H the left figure presents the 
initial configuration (in red) and the first iteration (in blue), the second shows 
the entire evolution for 100 iterations (unless stated otherwise), the last figure 
displays the scaled up configuration for the last few iterations. 
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Figure 4: A = i 




Figure 5: A = — i 




Figure 6: A = exp(i7r/4)/2 
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Figure 7: A = exp(i7r/4 + nr)/2 



4.2 Centroid gathering evolution and extensions 



As a second example, suppose that agent Vk is moving according to the following 
linear combination of its own position, the positions of agents higher in the hierarchy 
i.e. {Vk+i, ■ ■ ■ ,Vn-i} anci the positions of those lower than itself {Vq,Vi, ■ ■ ■ , Vk-i}'- 



N-l 



fc-1 



V k (t + 1) = aV k (t) + Pf Y W) + Pb Y W 



=fc+i 



or 



P(t+1) 



a ftp Pf ■ ■ ■ Pf 
Pb a p F ••• Pf 



P(i). 



Pb Pb a 

Note that if p F = Pb = (1 - at)/(N - 1), we will have 

1 - N 

V k {t + l) = aP k {t) + —^- Y 

JVa-1 / Na - I * 

-7\ i +1 — ' centroid 



N-l 



a — 1 



hence all agents move towards the time-invariant centroid on straight lines. 

For general Pf and Pb, the above matrix is /?s//?p-factor circulant and is diago- 
nalized by 

" 1 

(Pb/Pf) o 



P(i) = [FT] 







P(t), 



the modes or eigenvalues being given by 

N-l 



(j,i = a+y]p F 
k=l 



B 



fL 

Pf 



e -i w ki^ ; = ,...,AT-1. 
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Let us consider first the case of perfectly cyclic interaction, i.e., when (3g = ftp. In 
this case, the interaction matrix is circulant, and we have 



N-l 



= a + PFe~ i ^ kl , I = 0, . . . , N - 1 



k=i 



and 



Mo 
Mi 



a + (N- l)(3 F 



N-l 



a 



Pf+Y1 foe-^ 1 = a-0 F . 



k=0 



For normalization, we shall take (3f = (1 — o)/(N — 1) and then 

Mo = 1 

m = (Na-1)/(N -1), for all I. 



We now have that 



evolves according to 



P(t) = [FT]*P(i) 



P(t) 



t— »oo 



Hence 



P(t) 



jVa-1 
AT-1 



t— >oo 



[FTl 



1 




P(0) 



1 





[1,0,...,0]P(0). 



[1,0,...,0][FT]*P(0) 



N 



[1,1,...,1]P(0) 



i.e., as we have already seen, all points converge towards the centroid of the initial 
constellation. The convergence will be as follows: 



P^Wi-oo = P(t) 



[1,0,...,0]P(0) 



Na-1 
N-l 



P(0). 
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Figure 8: A = 1, a = 0.1, 100 iterations 



Therefore 



P"(t)^oo 



[FT1 



1 




[1,0, 



■0] 



[FT]*P(0) 



Not — 1 
N — 1 



P(0) 



V 



\ 



[1,1,...,1]P(0) 



/ 



Consequently, all agents will gather towards the centroid by moving on a line from 
P fc (0) to W) (see Figure EJ). 

Next suppose we have Pb ^ Pf- Then we have a A = Pb/Pf factor circulant and 
the modes of the P(i) evolution is controlled by 



in 



fe/AT 



e -^ 5 Z = 0,...,JV-1. 



Here 



Mo 



o 



a — Pf + Pf 



k=0 

Pb/Pf - 1 

{Pb/Pf) 1 '* - 1 



a: 



1 — a 
iV- 1 



+ 



1 — a 



A - 1 



iV - 1 7 V aV^v _ 1 



Similarly we have that 
W = 



a 



fc/n 



e 1 jv Kl 



1 — a 



a 



+ 



1 — a 



iV - 1 V N - 1 7 V M N e M l N -17' 



Ae 



-iff/ 



1 
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0.2 0.4 0.6 0.8 



Figure 9: A = 0.5, a — 0.1, 1000 iterations 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 





Figure 10: A = —1, a = 0.1, 100 iterations 



In this example too, as before, we have 
" 1 



P(t) 



t— »oo 



, N-l 

A" 



[FT! 



Mo 





and if is the dominant eigenvalue, we shall have 



[FTl 



x -i/n 



-(JV-l) 

A^iv-^ 



P(0) 



P(t) 



t^oo — Ho 



1 

N 



1 

X N-1/N 



[i,x~y N ,...,x- {N - 1)/N Mo) 



and depending on the values selected for A, we can get a wealth of interesting behav- 
iors while the solutions converge or diverge to infinity, displaying spiralling or in line 
marching. See Figures (H [TUl HH CGI [El [14] where we present a few interesting cases. 
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-0.005 0.005 0.01 0.015 0.02 0.025 0.03 



Figure 11: A = i, a = 0.1, 100 iterations 




Figure 12: A = — z, a = 0.1, 100 iterations 
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Figure 14: A = exp(z7r/4 + wr)/2, a = 0.1, 100 iterations 

5 Concluding Remarks 

We discussed in this paper a special type of cyclic multiagent interaction modeled 
by A-factor cyclic matrices. Such matrices allow explicit closed form diagonalizations 
via generalized Fourier transforms hence enable the analysis of the evolution of the 
swarm via a nice, geometric, modal decomposition process. It is expected that a 
wealth of further similar, structured and nearly cyclic interactions will also yield explicit 
closed form solutions for their asymptotic behavior. In fact, we may use evolutions 
that fix one, two [16] or several agents in the swarm and use circulant or A-circulant 
interactions for the rest of them leading to further highly structured matrices that can 
be diagonalized, and correspondingly leading to interesting and explicitly predictable 
and designable swarm dynamics. In closing, we note that Turing's morphogenesis may 
be regarded as a further example of such dynamics for points in the plane where the 
x and the y coordinates are subjected to different linear circulant transformations also 
readily generalizable to A-circulant maps [15] . An analysis of such swarm interaction 
for multiagent system is forthcoming. 
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